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ABSTRACT 

Many fast supernova remnant shocks show spectra dominated by Balmer lines. The Ha pro- 
hies have a narrow component explained by direct excitations and a thermally Doppler broadened 
component due to atoms that undergo charge exchange in the post-shock region. However, the 
standard model does not take into account the cosmic-ray shock precursor, which compresses 
and accelerates plasma ahead of the shock. In strong precursors with sufficiently high densities, 
the processes of charge exchange, excitation and ionization will affect the widths of both narrow 
and broad line components. Moreover, the difference in velocity between the neutrals and the 
precursor plasma gives rise to frictional heating due to charge exchange and ionization in the 
precursor. In extreme cases, all neutrals can be ionized by the precursor. 

In this paper we compute the ion and electron heating for a wide range of shock parameters, 
along with the velocity distribution of the neutrals that reach the shock. Our calculations pre- 
dict very large narrow component widths for some shocks with efficient acceleration, along with 
changes in the broad- to-narrow intensity ratio used as a diagnostic for the electron-ion temper- 
ature ratio. Balmer lines may therefore provide a unique diagnostic of precursor properties. We 
show that heating by neutrals in the precursor can account for the observed Ha narrow compo- 
nent widths, and that the acceleration efficiency is modest in most Balmer line shocks observed 
thus far. 

Subject headings: shock waves — acceleration of particles — ISM: supernova remnants 



1. Introduction 

Cosmic rays are widely believed to originate in 
supernova remnant (SNR) shock waves, because 
the cosmic-ray energy spectrum agrees with model 
predictions, because power-law distributions of en- 
ergetic electrons are seen in SNRs, and because the 
power required to maintain the cosmic ray popu- 
lation could be supplied by about 10% of kinetic 
energy of Galactic supernovae. The standard the- 
ory for the process is Diffusive Shock Acceleration 
(DSA), which is a first order Fermi process requir- 
ing that particles scatter between a gasdynamic 
subshock and plasma turbulence in a shock pre- 



cursor. Evidence for non-linear DSA comes from 



curved synchrotron spectra (IRevno lds & El lison 



19921 Kink et al. 1 120061 : I Allen et al. I l2008h 



dence for high compression factors ([Warren et al. 



evi- 



20051 : iCassam-Chenai' et al. 1 120081) and evidence 
for lower than expect e d downstream temp eratures 



(|Hughes et al. Il2000l: iHelder et al. 1120091) . How 
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ever, all this evidence is based on observation of 
downstream properties. The effects of precursor 
physics on the Ha emission described here offer a 
direct probe of the properties of the precursor. 

A crucial parameter for these models is the dif- 
fusion coefficient k, as it determines the precursor 
scale length, which is typically k divided by the 
shock speed Vs. Gas is compressed in the pre- 
cursor and accelerated to a fraction of the shock 
speed, and this compression is related to Vs, to 
the efficiency of particle acceleration and t o the es- 
cape of energetic particles from the region (|Bvkov 
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2005 ; Vink et al. 2010l) . Neutrals can impede 
the acceleration process by damping the turbu- 
lence nee ded to scatter particl es back to the shock. 
However, iDrurv et al. T (ll996l) found that the ac- 
celeration efficiency can be high as long as the den- 
sity and neutral fraction are not too large, though 
the maximum particle energy is reduced. 

One set of diagnostics for the physics of col- 
lisionless shocks is based on the emission from 
particles in the narrow i onization zone just be- 



hind a nonradiative shock ([Raymond II 199 It iHeng 
2010h . In particular, Ha photons from a non- 
radiative shock in partly neutral gas originate 
very close to the shock, and Coulomb collisions 
do not have time to erase such signatures as 
unequal electron and ion temper atures or non- 



Maxwellian velocity distr i butio n s (lLa ming et al 



1996 
2008 
seen 



Ghavamian et al 
2010h 



2001 



i s (E 
iRa; 



Raymond et al 
In the optical these shocks are 
Balmer line filaments whose pro- 



collisions with protons and electrons. If the den- 
sity is very high, neutrals and protons are tightly 
coupled by charge transfer. In that case, the 
neutrals are compressed along with the protons 
and adiabatically heated. They also share in any 
other heating of the protons, such as dissipation 
of Alfven waves generated by cosmic-ray stream- 
ing. On the other hand, if the density is very 
low, neutrals pass through the precursor and the 
shock without interacting at all, preserving their 
pre-shock velocity distribution. 

The intermediate case is more complex. A 
shock that efficiently accelerates cosmic rays is 
strongly modified, and gas reaches a significant 
fraction of the shock speed in the precu rsor 
( Vladimirov et al. l l2008HWagner et al. Il2009h . If 
neutrals and ions are fairly well coupled, they can 
be described as fluids whose relative speed gives 
a frictional heating sim ilar to that in C-shocks 



as pure 

files show a narrow component characteristic of 
the pre-shock kinetic temperature and a broad 
component closely related to the post-shock pro- 



( Draine fc McKee I I1993T ). If a neutral encoun- 



ton temperature ( Chevalier & Ray mond Il978t and become a pickup ion ([Raymond et al 



Heng||20ld: Ivan Adelsberg et al. 1120081) . The i 



ters this high speed compressed plasma without 
having been brought gradually up to speed by 
many previous charge tra nsfers, it can be ionized 

~ 20081: 



tensity ratio of the broad and narrow components 
is determined by the e lectron to ion tempera- 



Ohira fc Takahara 1 12010|) like those observed in 
the solar wind ( Moebius et al. Ill985l ). It can then 



ture ratio at the shock (Ghavamian et al. 2001 



van Adelsberg et al. 1 120081 iHelder et al. Il2010l ) 

The Balmer line profiles also contain signatures 
of shock precursors. In general, the narrow com- 
ponent line widths are 40 to 50 km s _1 , indicat- 
ing temperatures around 40,000 K. If that were 
the ambient ISM temperature, there would be no 
neutrals to create the Balmer line filament, so the 
width is interpreted as an indication of heating in a 
narrow pr ecursor too thin to completely ionize the 



have an energy on the order of 1 keV, which it can 
share with the other protons. Electron heating 
is more uncertain, b ut it can occur by m eans of 
Lower Hybrid waves ([Cairns fc Zank 1120021 ) . If the 
electrons are heated they can excite and ionize H 
atoms, changing the Ha profile and the broad-to- 
narrow line ratio used as an electron temperature 



diagnostic ([Ghavamian et al. 1 12001) 



hydrogen jSmith et al.lll994t iHester et al. Ifl99"3 



Lee et al. Il2007t ISollermanet al. 1120031 ). Faint 



In this paper we compute the proton, neutral 
and electron temperatures in the precursors for a 
variety of parameters, along with the ionization 
and excitation of H atoms. We consider the ef- 
fects of these process es on Balmer line d i agnos - 



emission ahead of the sharp filament is interpreted 
as emission from the compresse d and heated pre 



tics currently in use. lOhira fc Takahara I (|2010h 



curso r gas ([Hester et al. 
20101) . 



1994; iLee et al. I L2007 



This paper considers the role of neutrals in 
heating the precursor plasma and computes the 
properties of precursor Ha emission. While cos- 
mic ray pressure in the precursor can compress, 
heat and accelerate ions and electrons by means 
of plasma turbulence and magnetic fields, the neu- 
trals only interact with the precursor by means of 



considered the effects of neutrals on the velocity 
structure of the precursor, the compression ratio 
and the acceleration process. They found that 
the pickup ions can reduce the compression by 
the subshock and enhanc e proton injecti o n into 



Morlino et al~l (|2010l ) 



the acceleration process, 
self-consistently computed the particle accelera- 
tion and heating due to neutrals, but within the 
fluid approximation for both neutrals and ions. In 
this paper we emphasize the effects on the Ha line 
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profile. 

2. Model Calculations 

We parameterize the precursor structure in a 
relatively simple manner. We assume that the pre- 
cursor accelerates and compresses the interstellar 
gas over a length scale k/Vs, where n is the dif- 
fusion coefficient for cosmic rays near the cutoff. 
Effective cosmic-ray acceleration requires k on the 
order of 10 24 cm 2 s _1 , and estimates based on the 
scales of Ha precursors are 2 to 4xl0 24 cm 2 s _1 
(|Lee et al. 1120071 l2010h . We do not consider the 
second order effects of momentum and energy de- 
position by the neutrals on the precursor length 
scale. 

We assume an exponential form, so that the 
compression is given by 



X 



i + Oa-i) 



3 (xVs/k) 



(1) 



where x is negative ahead of the shock and xi is 
the compression ratio just upstream of the sub- 
shock. It is related to the fractional pressure 
of cosmic rays behind t he shock, w=Pcr/ (Pq + 
Pen), by equation 9 of IVink et al~l (|2010h . We 



simplify this equation with the assumption that 
for w<0.8 the compression in the gas subshock 
equals 4, so that 



Xi = (l-w/4)/(l-w) 



(2) 



Mass conservation implies velocities V — Vs fx m 
the frame of the shock. 

To compute the proton and electron tempera- 
tures we include adiabatic compression, Coulomb 
energy transfer between protons and electrons, en- 
ergy losses due to ionization and excitation of Hy- 
drogen, and heating terms. 

We assume that any neutral that interacts with 
the plasma at position Xi joins the proton flow at 
that position. If the interaction was charge trans- 
fer, a new neutral is formed with the bulk speed 
and thermal speed of the protons at xi. Thus the 
neutrals arriving at Xi are those that last went 
through charge transfer at all upstream positions 
Xj, and they have the speeds, Vj, of the plasma 
at Xj. Each ionization of a neutral from Xj at Xi 
deposits energy 0.5m p (vj — Vi) 2 . We assume that 
the energy is quickly thermalized amon g the pro- 
tons, unlike lOhira fc Takahara who as- 



sumed a pickup ion velocity distribution. The 
thermalization time scale is very uncertain, be- 
cause full kinetic calculations have not been car- 
ried out. However, in the highly turbulent pre- 
cursor there are many wave modes besides Alfven 
waves that can thermalize the protons, in partic- 
ular those associated with bump-on-tail, mirror 
and firehose-like instabilities ( Winske et al. Ifl985t 



iGarv Ill978t ISagdeev et al. Ill986l ) We also 



ignore 

heating due to Alfven wave damping or shocks ex- 
cited by the cosmic-ray pressure gradient in order 
to isolate the effect of the neutrals. Therefore, we 
compute a lower limit to the heating. 
For electrons, we follow Cairns fc Zank 



(2002), 



who found that ionization of fast neutrals forms a 
ring beam, in which all the particles gyrate around 
the magnetic field with the same speed but differ- 
ent phases. The ring beam is unstable, and pro- 
vided that the beam velocity (in this case the rela- 
tive velocity of bins i and j) is less than 5 times the 
Alfven speed, it transfers a significant amount of 
hea t to electrons via Lower Hybrid waves. We fol- 
low Cairns fc Zank ( 20021 ) in taking this fraction 
to be 10%. Again, to isolate the effects of neu- 
trals we ignore any heating of electrons by Lower 
Hybrid waves generated by cosmic-ray streamin g 



(|Ghavamian et al. |[2007tlRakowski et al. Il2008h . 



Charge transfer rates are taken from (jSchultz et al 



2008) using the quadrature sum of the thermal 
speed and the ion-neutral relative speed. Ion- 
ization and excitat ion rates are compute d from 
cross sections from Janev fc Smith ] (<1993l) by in- 
tegrating over the electron velocity distribution 
including the relative electron-neutral flow speed. 

3. Results 

Figure 1 shows a set of models for a shock speed 
of 2000 km s" 1 with k = 2.0 x 10 24 cm 2 s" 1 , a pre- 
shock density of 0.2 cm -3 and a neutral fraction 
of 0.2. The four models have ratios of cosmic-ray 
partial pressure to total pressure behind the sub- 
shock, w, of 0.1, 0.3, 0.5 and 0.7. The compres- 
sion ratios just ahead of the subshock are 1.0833, 
1.3214, 1.75 and 2.75. In the high V s , high ef- 
ficiency models the neutrals are not compressed 
to this level, because of collisional ionization and 
because some pass through without charge trans- 
fer. The protons and electrons are strongly heated 
in the more efficient models, but the electrons are 
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much cooler than the protons. The drop in heating 
just before the subshock in the 70% efficient model 
results from the reduced number of neutrals. 

Figure 2 shows the velocity distributions of the 
neutrals perpendicular to the shock just before the 
subshock. Note that the w=50% model shows a 
narrow component due to neutrals that last ex- 
perience charge transfer far upstream, along with 
a broader component of particles that undergo 
charge transfer close to the subshock. 

Figure 3 shows a grid of models of the neu- 
tral velocity distribution at the shock for a range 
of shock speeds and cosmic-ray partial pressures. 
The panels show the FWHM measured directly 
from the computed velocity distribution and the 
kurtosis, which would be 3.0 for a Gaussian distri- 
bution. Kurtosis is a problematic statistical mo- 
ment for real data because it is sensitive to noise 
far from the line center and the choice of back- 
ground level. However, for the theoretical profiles 
computed here it highlights cases in which some 
neutrals undergo charge transfer close to the sub- 
shock and others do not. We also show the frac- 
tion of incident neutrals that survive up to the 
subshock and the average number of excitations 
to the n=3 level per incident H atom. Not all of 
these excitations will result in Ha photons because 
some Ly/3 photons escape, but this is a convenient 
comparison to the 0.2 to 0.25 Ha photons per H 
atom produced in the post-shock region. 

4. Discussion 

The interaction of neutral hydrogen with the 
ionized plasma in cosmic-ray precursors described 
above offers an important tool to measure the 
properties of cosmic-ray precursors. The out- 
come of DSA is very much influenced by phys- 
ical processes in the precursor, which are not 
well determined. For example, non-adiabatic 
heating and magnetic field amplification due 
the presence of cosmic rays tend to decrease 
the overall compression facto r from X12 >> 20 



energy escape automatically implies a lower do wn- 



(e.g 
to 7 



Berezhko fe Ellison 



< 



J02. 



ICaprioli et al. 



< 1 5 
2008). 



1999; Blasi et al. 



( Vladimirov et al. 



2005) 



2008; 



In addition, if the Alfven 
waves in the precursor have some drift velocity 
this will affect the cosmic-ra y pressure profile 
( Zirakashvili fc Ptuskin 20081) . which limits the 
escape of energy from the shock region. A lower 



stream cosmic-ray pressure rtVink et al. Il2010h . 

As shown here, neutr als will influen c e the 
physics of the precursor. Morlino et al. ( 2010l ) 



treated the neutrals as a fluid coupled to the ions 
by charge transfer for a unified treatment of the 
heating and dynamics of the precursor, but the 
fluid approximation is only appropriate if neutrals 
and ions are coupled fairly well. They obtained 
a FWHM of 46 km s" 1 for the Ha line in a 2000 
km s _1 shock with modest efficiency, a pre-shock 
density of 1 cm -3 and 50% neutral fraction. For 
similar parameters we find a non-Maxwellian pro- 
file with smaller FWHM and broader wings. 

Neutrals can also damp plasma waves, which 
limits the efficiency of cosmic-ray acceleration. 
This damping is caused by the central processes 
described above: ch arge exchange and ionization. 
Drurv et al. ( 19961 ) found that the maximum par- 



ticle energy, and therefore the maximum acceler- 
ation effic iency, is considerably hig her than sug- 
gested by iDraine fc McKee" ( 1993 ). In addition, 
the heating due to neutrals penetrating the pre- 
cursor is a form of non-adiabatic heating. Energy 
dissipated in the precursor limits the amount of 
free energy available for shock acceleration. If 
the neutrals ionized in the precursor behave as 
pickup ions rather than thermalizing with the pro- 
tons, the injecti on efficiency and particle s pectrum 



will be affected (jOhira fc Takahara" 2010l ) . In any 



case, the heating of electrons in the precursor is 
poorly know, and that will stongly affect the ion- 
ization and excitation of H atoms, which in turn 
will affect the intensity ratio of the broad and nar- 
row components as well as the narrow component 
line width. 

The physics of neutral-ion coupling means these 
processes are not only sensitive to cosmic-ray pres- 
sure and the structure of the precursor, but also 
to the pre-shock density and neutral fraction. For 
example, the protons and neutrals in the Cygnus 
Loop nonradiative shocks ( Salvesen et al. I 20091) 
are tightly coupled, and they behave nearly adi- 
abatically, while the pre-shoc k density in SN1006 
is so low ( Acero et al. 1 120071) that neutrals pass 
straight through it. Th is may explain the nar - 



row line width seen by ISollerman et al. ( 2003 ) 



in SN 1006 in comparison to broader narrow lines 
observed for other young SNRs. The number of 
charge transfer events for an average neutral in 
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the precursor can be estimated roughly as 



N C t = nLaxi/Vs, 



(3) 



where L is the precursor length scale, and the 
charge transfer cross section, a, declines slowly 
with velocity below about 2000 km s _1 , then very 
rapidly. 

Comparing our results to observations, it is ob- 
vious that the observed narrow-line Ha widths are 
in general smaller than predicted by our calcula- 
tions for efficient shock acceleration (i.e., w > 0.5). 
This may indicate that none of the shocks inves- 
tigated so far accelerate particles efficiently. How- 
ever, more work is needed before such a conclu- 
sion can be drawn, as the line width depends 
also on pre-shock density and shock velocity. For 
very high shock velocities combined with low den- 
sities the neutrals hardly interact in the precur- 
sor, leading to narrow line widths. Such may 
be the case for the northeaste rn re gion of RCW 

(|2009li reported a 



86, for which iHelder et al 



high cos mic-ray acceleratio n efficiency (w > 0.5 



see also Vink et al. 2010( ). For this region the 



pre-shock density may be as low as n < 0.1 
( Vink et al. ] l2006h . which, combined with the 
high velocity (V s > 3000 km/s), gives few inter- 
actions in the precursor and widths < 100 km/s. 
Note that this is smaller than could be measured 
given the moderate spectral resolution of the mea- 
surement. 

Perhaps the most striking result from these 
calculations is that for efficient shocks near 1000 
km s _1 a substantial number of Ha photons will 
be produced in the precursor. These will usu- 
ally be included in the narrow component, po- 
tentially affecting the electron temperature es- 
timate base^_oji_the__broad^to-narrow intensity 



ratio (jGhavamian et al. I 120011 ). Narrow com- 



ponent emission from the precursor could ex- 
plain the broad-to-narrow intensity ratios that 
cann ot be fit by models of post-shock emis 



sion ( van Adelsberg et al. 20081 : Rakowski et al. 
2009). In extreme cases, emission from the precur 
sor might also contribute to the broad component, 
possibly accounting for t he non-Maxwellian pro - 



file seen in Tycho's SNR (jRavmond et al. I l2010h 



and generally leading to an underestimate of the 
shock speed. Both of these conclusions depend on 
the diffusion coefficient and the electron heating, 



however. 

Other Ha line measurements show that the nar- 
row lines are broader than one might expect for 
temperatures of typical HIT regions, but smaller 
than 50 km/s (jSollerman et al. H2003h . (Not all of 
the narrow line emission comes from the precur- 
sor, but the narrow line emission downstream is 
determined by the velocity distribution in the pre- 
cursor.) Another effect of charge exchange in the 
precursor is that neutrals enter the downstream 
shock region with a velocity offset with respect 
to the local interstellar medium, as seen in Ty- 
cho's supernova remnant (Lee et al., 2007). For 
shocks observed face on this should produce a nar- 
row line offset, which for the combined front and 
back side of the remnant should lead to two nar- 
row lines. The spec tra of several LMC remnants 
(ISmith et al. Ill994h do not show such an effect. 
For one of the remnants in this set, SNR 0509- 
67.5, the cosmic-ray ac celeration efficiency was es- 



timated to be w « 0.2 (jHelder et al. Il2010t ). 



We note that several improvements should be 
made to the calculations presented here. Addi- 
tional heating due to wave dissipation can heat 
the protons, resulting in larger narrow component 
line widths, or it can heat electrons, increasing 
the Ha narrow component intensity and reduc- 
ing the number of neutrals that reach the shock, 
especially if the electron velocity distr ibution is 
non-Maxwellian ( Laming fc Lepri"ll2007l) . In addi- 
tion, ionization and excitation by proton and he- 
lium ion impact are importan t at high relative ve- 



ipc 

locities ( Laming et al. I 19961 ). and at high shock 



speeds the velocity distributions of particles are 



ot pa 

anisotropic dHeng fc McCrav 1 120071 : LHeng et al 



20071 : Ivan Adelsberg et al. 1120081) . Amplification 



of the magnetic field may also be important, and 
radiative transfer calculations in the Ly/3 line must 
be done to compute the Ha emission. We plan to 
address these issues in future work. 
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the VIDI grant awarded to JV by the Netherlands 
Science Foundation (NWO). 
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Fig. 1. — Plots of a) proton temperature, b) electron temperature, c) neutral density and d) heating rate 
for for models having n = 2.0 x 10 24 cm 2 s _1 , pre-shock density = 0.2 cm~ 3 and neutral fraction 0.2. The 
models assume a post-shock cosmic-ray pressure of 10%, 30%, 50% and 70% of the total pressure, with the 
10% curves at the bottom and the 70% curves at the top in all four plots. 
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Fig. 2. — Velocity distributions perpendicular to the flow direction at the gasdynamic subshock for the four 
models shown in Figure 1. Note the narrow central component and broader wings in the model with 50% 
cosmic-ray pressure and the large width predicted by the 70% cosmic-ray pressure model. 
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Fig. 3. — Model grid for k — 2.0 x 10 24 , a pre-shock density of 0.2 cm -3 and neutral fraction of 0.2. Models 
are shown for Pcr/(Pcr + Pg) = 0.1 to 0.8 at the shock front. The panels show the FWHM and kurtosis 
of neutrals that reach the subshock, the fraction of neutrals that reach the subshock and the number of 
excitations to n=3 in the precursor per incident neutral hydrogen atom. As in Figure 1, the 80% models are 
the extreme cases. 
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